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Abstract 

In the collinear factorization of the form factor for the transition 7*71" vr the hard part 
contains double log terms. These terms will spoil the perturbative expansion of the hard part. 
A simple exponentiation for resummation leads to divergent results. We study the resummation 
of these double log's. We make an analysis to show the origin of the double log's. With the 
understanding of the origin one can introduce soft factors and Nonstandard Light-Cone Wave 
Functions(NLCWF) to derive a factorized form for the form factor, where the hard part does not 
contain the double log's. There is a perturbative relation between NLCWF and the standard 
Light-Cone Wave Function(LCWF). Beside the renormalization scale ^ the introduced NL- 
CWF's and soft factors have extra scales to characterize the double log's. Using the evolutions 
of the extra scales and the relation we can do the resummation of the double log's perturba- 
tively in sense that LCWF's are the only nonpertubative objects in the resumed formula. Our 
results with some models of LCWF show that there is a significant difference between numerical 
predictions with or without the resummation. 
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1. Introduction 



Although the interaction of QCD is asymptoticahy weak , the perturbative theory of QCD can 
not directly be used to study hadronic processes involving large momentum transfers because of 
quark confinement. One needs to separate or factorize long-distance- and short-distance effects. 
Only the latter, which are characterized by large momentum transfers denoted generically as Q can 
be studied with perturbative QCD. For an exclusive process it has been proposed long time ago 
that such a process can be studied by an expansion of the amplitude in 1/Q, corresponding to an 
expansion of QCD operators in twist [H [2]. The leading term can be factorized as a convolution of 
a perturbative coefficient function and light-cone wave functions of hadrons. The light-cone wave 
functions are defined with QCD operators. The perturbative coefficient function, which can be 
safely calculated with perturbative QCD, describes hard scattering of partons at short distances. 
This is the so-called collinear factorization. In this factorization the transverse momenta of partons 
in parent hadrons are also expanded in the hard scattering part and they are neglected at leading 
twist. 

Perturbative coefficient functions in collinear factorization contain large logarithms at higher 
orders of a^. These large logarithms are dangerous and can spoil the expansion in the sense that the 
expansion does not converge. A resummation of large logarithms is often needed to have a reliable 
prediction. In this paper we study the resummation in the process 7*7r+ — > tt^. The process is 
described by the electromagnetic form factor F[Q) with Q as the virtuality of the virtual photon. 
For large Q the form factor can be written as a convolution 



F{Q) ~ 0(x) (t>{y) ® H{x, y) 1 + O(AVQ^) 



with A as a scale characterizing nonperturbative effects. In the above (t){x) is the light-cone wave 
function(LCWF) of vr with the momentum fraction x carried by a parton, H is the hard part which 
can be calculated safely as an expansion in q^. The one-loop correction to the hard part has been 
studied in [3llll[5l[6l[7]. It has been found that the hard part at one-loop contains double log terms 
like In^ x and a1 \t? y. At n-loop level the hard part will contain aj^" In^" x or a^"^" In'^" y. 
Those double log terms will become divergent when x or y is approaching to zero and can spoil 
the perturbative expansion of H. Although the double log terms are integrable with x or y in the 
convolution and give finite contributions, but they are significant corrections. The purpose of our 
study is to resum those double log terms. 

We study the resummation of the above double log terms in the collinear factorization by re- 
factorizing the double log terms in the hard part. For doing this we have to understand the origin 
of the double log. Given the factorized form in the above, H will receive in general contributions 
from the form factor and the LCWF. If we use a finite quark mass to regularize the collinear 
singularities, we can show that a part of double log terms come from LCWF. This part can be 
re-factorized by using a nonstandard light-cone wave function(NLCWF) which has a perturbative 
relation to LCWF. The remaining double log terms come from the form factor. They can be re- 
factorized into soft factors. In introducing NLCWF and soft factors non-light like gauge links play 
an important role. The importance of using non-light like gauge links has been shown in Transverse- 
Momentum-Dependent(TMD) factorization for inclusive processes [8l [9l [TOl 111]. With these gauge 
links extra scales are introduced to control these double log terms. Using evolution equations of 
these scales we are able to resum the above double log terms. The outlined approach has been used 
to resum double log terms in 7r7* — > 7 in [12j, where the resummation has a significant effect. Our 
approach is similar to the threshold resummation in inclusive processes studied in |13] . 
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The resummation of these double log terms has been studied in [11] with the /c-r-factorization. 
The /cy-factorization has been widely used in studies of i?-meson decays. However, such a fac- 
torization is not gauge-invariant because hard parts are extracted from scattering amplitudes of 
off-shell partons. The amplitudes of off-shell partons are in general gauge-dependent. Recently it 
has been shown that the hard parts in the factorization receive at loop-level divergent contribu- 
tions which are gauge-dependent [15j ■ One may find further discussions of the issue in [IB]. The 
A:T-factorization takes transverse momenta of partons into account. These momenta are neglected 
at leading twist in the collinear factorization. It should be noted that the effects of the transverse 
momenta can be taken into account in a gauge-invariant way by using the so-called Transverse 
Momentum Dependent (TMD) factorization [T71 H]. It is possible to use the TMD factorization to 
perform the resummation. But it is at moment not possible to give numerical predictions in detail 
because the involved nonperturbative objects are poorly known. 

Our paper is organized as the following: In Sec. 2 we introduce our notations and give a brief 
discussion about the consequence of the double log terms from higher orders. In Sec. 3 we explain 
the origin of double log terms in the hard part. In Sec. 4 we introduce NLCWF and soft-factors to 
re-factorize double log terms. In Sec. 5 we show that these double log terms can be resummed and 
give our numerical results. Sec. 6 is our conclusion. An appendix is given to discuss the problem 
of gauge invariance in A;r-factorization for the case studied here. 

2. Notations 

The electromagnetic form factor of 7r+ is defined as: 



We will use the light-cone coordinate system, in which a vector is expressed as = (a^, a ,a±) = 
((a° + a3)/V2, (a° - a^)/V2,a\a'^) and = (a^)'^ + (a^)^. Two vectors = (1,0,0,0) and 
= (0, 1, 0, 0) are introduced. We take a light-cone coordinate system in which the momenta are 



P'^ « (P+, 0,0,0), K^' ^ {0,K-,0,0), = K - P, = -q'^ ^ K~ . (3) 



When is very large, the form factor takes a factorized form. To derive the factorized form, i.e., 
to determine the hard part, one usually replaces hadronic state with partonic states to calculate the 
form factor and LCWF's. We replace the initial state |7r"'"(P)) with the partonic state \d(pi)u{p2)) , 



= (xoP+, 0,0,0), = (xoP+, 0,0,0), k>^ = {0,yoK-,0,0), k!^ = {0,yoK- ,0,0) (4) 



with xq = 1 — xq and yo = ^ ~ Uo- We take massless quarks and all singularities are regularized 
with dimensional regularization. 

The contributions to the form factor of the replaced partonic states with the virtual photon 
attached to the u-quark are from two diagrams given in Fig.l. To reduce the number of dia- 
grams which need to be calculated, we take /i = — . Fig. lb will not contribute for /i = — . It is 
straightforward to obtain the u-quark contribution: 



{7T+{K)\J^^\7r+{P)) = F^{Q){P + K) 



(2) 



given as: 



and the final state {■7t~^{K)\ with the partonic state {u{k2)d{ki)\. The quark pairs are in color- 
singlet. The momenta are: 



la 



u 




(5) 
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(a) (b) 

Figure 1: The leading order contributions from the u-quark. The black dot denotes the insertion 
of the electromagnetic current, i.e., inserting a 7^. 

where e„ is the electric charge of the u-quark. Similarly, one can obtain the contribution from the 
d-quark through the symmetry of charge conjugation. 
The definition of LCWF is: 

^{x,^i) = y'^e-^^^-(0|J(0)Lt(oo,0)7+75in(oo,z)tx(z-n)|^+(P)), (6) 

where q{x){q = d,u) is the light-quark field. L„ is the gauge link in the direction n: 

Ln{oo, z) = P exp ^—igs J dXn ■ G{Xn + z)^ . (7) 

If we replace vr with the parton state we have at the tree level for the wave function: 

(^(°)(x,;u) = 5(x-xo)0o + ---, 'Po = HPih''lMP2)/P^, (8) 

where • • • stand for the states of the quark pair with quantum numbers other than that of vr. For 
the outgoing ir^ LCWF is defined with quark fields which are separated in the direction I and we 
the corresponding gauge link is along the direction of I. With the definitions one easily obtains the 
f actor ized form: 

SttOs 1 

-^y 
n{x,y,Q,fi) = l + 0(a,). (9) 

In the above the hard part TC can be calculated as an expansion in as- The correction to the 
factorized form factor is power-suppressed. 

At one-loop level the hard part receives corrections which have double log terms [3l HI EJ [6l [7] : 

W(x,y,g,/x) = 1 + ^^ [in^x + ln^yj +..., (10) 

and at higher orders of a" those terms like a" In^" x and a" In^" y, which can lead to an divergent 
series in a^. If we take asymptotic form of LCWF, i.e., (j){x) ~ x(l — x), one easily finds that the 
1 + n-th order contribution to the form factor: 

~ a,^+"(2n)!. (11) 

The factorial increase of the 1 -|- n-th order contribution can spoil the perturbative expansion. 
Therefore, these double log terms need to be resummed in order to have a reliable expansion. It 
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should be noted here that a simple exponentiation of these double log's leads to a divergent form 
factor because the coefficient in the front of the double log's is positive. 

For the fcr-factorization in the case studied here the same diagrams in Fig.l gives the leading- 
order result, in which the partons are off-shell. Therefore, the result is not gauge-invariant. We 
will discuss this in detail in Appendix. 



3. The Origin of the Double Log Terms 



In this section we analyze the origin of the double log terms in the hard part. The one-loop 
correction has been calculated in[3l[4l[5l[6l[7j. In[3l[7j the individual contributions from each 
one-loop diagram are listed in detail. From their results we find with /i = — in Eq. (2) that the 
diagrams given in Fig. 2 give contributions containing double log's. In order to have short notations, 
we denote in this section the momentum fraction xq and yo in Eq.(4) as x and y, respectively. 

In general the double log terms are generated if the integration region of loop momenta has an 
overlap between collinear- and infrared region. The infrared region always exists. The collinear 
region can exist in two cases. One case is that a loop momentum k is directly collinear to an external 
on-shell particle. In this case the integration of k will generate a collinear singularity regularized 
with a mass of the external particle or with dimensional regularization. Another case is that in 
some kinematical region, e.g. when x or y approaches zero, some propagators can become nearly 
on-shell and the collinear region is formed when k is collinear to the momenta of these propagators. 
The integration over k will generate the corresponding collinear singularities appearing in the form 
as Inx or Iny. If there is an overlap of these collinear regions with the infrared region, double log 
terms will be generated. We will call those double log terms in the former case as Type I double 
log's, and these in the later case as Type II double log's. The existence of the two types of double 
log's will be explained in detail with the examples in Fig. 2. 

Let us first consider the contribution from Fig. 2a. This contribution contains double- log terms 
like In^ x and In^ y. It also contains a double pole term. These terms come from the overlap region 
of collinear- and soft gluons. The double-pole and double-log terms can easily be found by the 
approach of soft gluons where one neglects in the first step all loop momenta k in the nominator 
and picks the dominant terms in the denominator for — > 0. Then we have the soft part from 
Fig.2a: 



FniQ) 



2a " \2Nj J (27r)4((pi + A:)2 + ie)((A;i+fc)2 + ie)(A:2 + ie)' 




Api ■ ki 




(12) 



The reason for the double log's and the double pole in Fig. 2a is the following: The exchanged 
gluon 1 can be collinear to the initial q with the momentum pi, it can also be soft. This overlapped 
region produces an double pole and the double log of pf , i.e., In^ x. Similarly, the gluon can also be 
collinear to the antiquark in the final state with the momentum ki. Again there is an overlapped 
region of the soft momentum. This produces another double pole and the double log of k^ , i.e., 
In^y. The sum of the two double pole contributions is given in the above. It is clearly that this type 
of double log's is associated with partons in the initial- or final state. According to the discussion 
at the beginning of this section, these double log's are Type I double log's. Actually, the double 
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P2 



Pi 



(a) 







(/) 



Figure 2: One-loop correction with double log terms, the original labels of the diagrams in [7] are 
also given. 



log's and the double pole contributions from Fig. 2a are canceled in the final result because the 
quark pair in the initial- and final state is neutral in color. We note that the soft gluon emitted 
by the initial q can be absorbed by the final q and the final q quark, the former contribution is 
given by Fig. 2a, the later contribution is given by Fig.2f when the gluon 2 is soft. The sum of the 
two contributions represents the interaction of the soft gluon with the quark pair in the final state. 
The interaction of the soft gluon with the quark pair vanishes because the pair is color-neutral. 
Therefore, the double pole and double log associated with pi in Fig.2a are canceled by those in 
Fig.2f. With the same reason, one finds that the double pole and double log associated with ki in 
Fig. 2a are canceled by the soft contribution from Fig. 2d with the gluon 1 as the soft gluon. We 
will show this explicitly. 

We turn to the contribution from Fig. 2b. In Fig. 2b any one of the two gluons from the gluon 
splitting can be soft. We denote the momentum of the gluon 1 as k. The gluon 2 carries the 
momentum (ki — pi — k). If the gluon 1 is soft, we have the approximation for x ^ 0: 



„4 fabc 



(fk 



8/ci ki ■ p2 



K~{pi - kif{P - kiY J (27r)4 (fc2 + ie){{ki - pi - k^ + ie) 

■—^ \uik2)rT'T''u{p2Mpi)rT%{ki)] +■■■ 

— 2^2 «~ + L J 



1 



h 7 -|- in \ mx m x 

e 47r(u2 / 2 



(13) 



The term 2p^/c~ -|- ie) is the quark propagator with the soft gluon approximation. This 
contribution only generates the double log of x. It also contains a collinear singularity from the 
region in which k is collinear to p2- For the case that the gluon 2 is soft, no double log term is 
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generated. There is also no double log term In^ y for y ^ 0. For our purpose it is important to 
understand the above double log terms. 

In Eq.(13) there are two terms with the double log oi pf . One is associated with a collinear 
divergence, one is given as — (ln^x)/2. These two terms have different origins. The initial quark 
emits a collinear gluon and the gluon is absorbed by the virtual gluon attached to the antiquark 
line. The virtual gluon becomes on-shell when x approaches to 0. With x = one will find a double 
pole, one is collinear, another is I.R. one. In the case with x ^ 0, the I.R. singularity is regularized 
and it generates one \nx. Another Inx is generated from the collinear region. This is the origin 
for the first double log. It is clear that this double log is associated with the collinear region of p2 
and this double log is a Type I double log. The existence of Type I double log also depends on 
how the collinear singularity is regularized. To see this more clearly, we can use a quark mass m 
to regularize the collinear singularity. Then Eq.(13) becomes for x ^ 0: 



A fabc 



9tf 



8k^ ki ■ p2 



K~{pi - kiY{P - kiY J (27r)4 (fc2 + ie){{ki - pi - kf + ie) 

■- [u(A;2)7°T^T^(p2)^)(pi)7"r"r;(fci)l +•• 

\k — P2) — m + «e L J 



OsNc 

An 



1 "I , - 1,2- 

m —-77 m X m x + 

yQ^ 2 



+ ••• 



(14) 



Prom the above, one can see that the discussed double log disappear with the quark mass. It should 
be noted that the hard part does not depend how the collinear singularities are regularized. If one 
uses the mass to regularize the collinear singularity, one should also calculate the wave functions 
with the mass in order to subtract collinear singularities in the scattering amplitude. Then the 
Type I of double log's appearing in the hard part come from the wave functions. As we will show 
that Type I double log's can be factorized with NLCWF introduced in [T^ . 

The second double log in Eq.(13), or that in Eq.(14) has a different origin than that of Type I 
double log's. The virtual gluon attached with the antiquark line carries the momentum pi — ki. For 
X — > 0, the virtual gluon approaches on-shell, its momentum is nearly collinear to —ki. Here again 
the gluon 1 can have a momentum region which is the overlap region with the loop momentum is 
collinear to —ki and soft. For small non-vanishing x the collinear singularity and the I.R. singularity 
are simultaneously regularized by x, and it results in the second double log In^ x. This double log 
is a Type II double log. As we will show later, this type of double log's can be factorized by certain 
soft factors. 

The contribution from Fig.2e contains the double log In^ A;j~, i.e., In^ y. This double log can be 
obtained by making the approximation for the contribution in the limit of A; ~ 5 and fci ~ 5 with 
(5 — > 0, where we denote the momentum of the gluon 2 as k. Under the approximation we obtain: 



2e 



9 s 



f 



abc 



K- {pi - kiY 

d^k 2p+ 



u{k2h''T''T''u{p2)v{pih"T%{ki) 



(2^)4 P+ (ik2 + ie){{ki -pi- ky + ie){k+ + ie) 



OsNc 



In y h 7 — In 47r I + In y In 



+ 



+ In y ~ 2 ^ 



+ ••(15) 



Again, the first double log of y is Type I double log, and the term — (In^ y)/2 is Type II. Analyzing 
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the contribution from the region of pi ~ /c ~ 5 and the cases with the soft gluon as the gkion 1, we 
do not find any other double log's. 

A complicated case is the contribution from Fig.2f. We first consider the case where the gluon 
2 with the momentum k is soft. The contribution in the region of pi ~ A; ~ (5 we have the 
approximation : 



FAQ)\2f 



9s 



2NcK- 



u{k2)-i^T^u{p2)v{pi)-ipT^v{ki) 



2/un 



P-ki 



d'^k 



(27r)4 (fc2 + 



ie){{k + ie) 



+ ■ 



a. 




- +- 



"7 + In 



xyQ"^ 



2 xyQ^ 



+ 



(16) 



This contribution generates the double log of l'a^{pf), it also generates a double log of ln^(A;^). The 
later is irrelevant here for our purpose. From the above one can see that the double pole and the 
associated double log In^ x, i.e., In^(p^), is canceled by those from Fig.2a, as discussed before. The 
contribution from Fig.2f also contains a double log of k^ . To analyzing the contribution from the 
region of A; ~ A;i ~ 5 one can use the trick in [3j to decompose the production of five propagators 
into a sum of products of four propagators. Then one can find the dominant contribution from the 
region A; ~ fci ~ 5 as: 



FAQ)\ 



2/ 



NcK-pi ■ ki 

d^k 



u{k2)l^T^u{p2)v{j)i)^pT^v{ki) 

1 



(27r)4 (A:2 + ie){{k - k2f + ie)((P - ki - kf + ie) 



+ 



Ois 



+ In^ y - - In^ y 



+ 



(17) 



Again, the first double log of y is Type I double log, and the term — (ln^y)/2. For the case that 
the gluon 1 is soft with jU = — one does not find the double log of In^ A;j~. Hence, with = — only 
the soft gluon 2 generates the wanted double log's. The above results of Fig.2f are in agreement 
with those in [3l [7] . 

Performing a similar analysis of contribution from Fig. 2d, we find all relevant double log's from 
the region where the gluon 1 is soft. We denote k as the momentum of the gluon 1. In the region 
fei ~ A; ~ J we have the approximation: 



2d 



2x 



NcxQ^ 
d^k 



[u{k2)rT''u{p2)viPlhaT%iki)] 



(2^)4 (A:2 + ie){{k + kif + ie)({p2 + k^ + ie) 



Fi'\Q) 



a. 




- +- 



"7 + In 



+ iln^i4 
2 xyQ^ 



+ 



(18) 



Comparing the result in Eq.(12) for Fig. 2a one realizes that the double log In^y in the above is 
canceled by that in Eq.(12). The contributions from Fig. 2c does not have relevant double log's. 
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Adding the soft contributions from Fig. 2a, Fig. 2c, Fig. 2d and Fig.2f., all double poles are canceled. 
For the region pi A; ~ 5 we similarly have: 



2d 



NcQ^xyQ^ 



(27r)4 (A;2 + ie)((p2 + kf + ie){[pi - k^ - kf + ie) 



+ 



a. 



AnNr 



mx \ |-7 + ln 



^ 77 + In - - In X 
AiTfi^ I 2 , 



+ 



(19) 



Again, the first double log of x is Type I double log, and the term — (In^ x)/2. 

From the above analysis, one can find that the final contributions of double log's to the pertur- 
bative coefficient H only come from Type I and Type II double log's in Fig. 2b, Fig.2e, Fig.2f and 
Fig.2d. Adding every double log's together we have: 



n{x,y,Q,fi) 



^ I c^sjN^ - 1) 



(In^ X + In^ y) ~ 2 (^'^^ ~'~ ^^"^ 



-\ In x In y H 

Stt 47r 



R 1 ^ 



-(3 + Inx + Iny) In 



+ 



(20) 



where • • • stand for non-log, i.e., rational terms. In the first line the double log's given in [• • •] come 
from Fig. 2b, Fig.2e, Fig.2f and Fig. 2d. The Type I double log's are given in the first (•••), the 
Type II double log's are given in the second (•••). The results in the second line are from previous 
results in [Sj HI O O [7] . With the above results the double log In^ x from Fig.2b and Fig. 2d in 
Eq.(19) and the double log In^ y from Fig.2e and Fig.2f in Eq.(17) need to be resummed. 



4. Factorization of the Double Logs 

4.1 Factorization of Type I double log's 

We have seen in the last section that the origin of Type I double log's depends on how the 
collinear divergence is regularized. As discussed in detail in [12j, this type of double log's can be 
factorized by introducing NLCWF. The definition of NLCWF is: 



dz- (0|g(0)4(oo, 0)7+75^.(00, z-n)<7(z"n)K0(P)) 



2^ S{z-,Cu) 
Lu{oo,z) = Pexp(^-igsJ^ dXu ■ G{Xu + z)^ , tt'^ = (n+, 0, 0). (21) 

In the above the limit u~ ^ u~ should be taken. The main difference between LCWF and NLCWF 
is the gauge link in definitions. In LCWF the gauge link is obtained from La by setting = 0. 
It should be noted that the limit m"*" ^ or n~ ^ is taken after integrations of loop-momenta 
in the definition[8l [U \T0\ lllj . Otherwise the above definition is reduced to that of LCWF. The 
defined NLCWF depends on an extra parameter 

c?.™!.l(5^. (22) 
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The factor S{z , (^n) is the vacuum expectation value of a product of four gauge hnks. Details can 
be found in [12]. In [12] it has been shown that there is an interesting relation between LCWF and 
NLCWF. It reads: 



(a^,C,/^) 



dyC{x,y,C,n)(l){y,fi), 



(23) 



where the function C can be calculated perturbatively. The relation is in fact a factorization 
relation. NLCWF contains the same collinear singularities as LCWF. Therefore the function C 
does not contain any soft divergence. Detailed results of C at one-loop level can be found in |12j . 
For our purpose to do the resummation at leading order, it is enough to consider the function: 



y 



-C{y,xX,lj)- 



(24) 



From [12j we have: 



as{li){Nl 



+2^"C^ 



AttN^ 

2 1 



In^ X + In X In ^ + 2Li2 {x) 



vr 



vr 



+ - + 2\+0{al). 



(25) 



For vr in the final state moving with the momentum = (0, ii'~,0,0) one can introduce the 
NLCWF (j)(xXvifJ') with the replacement u ^ v with = {v~^ ,v~ ,0,0) and ^ v~ . The 
parameter Ci, is defined as: 



y2 



(26) 



Because NLCWF and LCWF contain the same collinear singularities, one can also use NLCWF 
to factorize the form factor. With these NLCWF's the form factor can be factorized as: 



Siras 



^^4){x, Cu,fx)4>{y, Ci},n)n{x, y, Cu, Cv, Q, At), 
X y 



Comparing two factorizations we have the relation: 

'H{x,y,Cu,Cv,Q,l^) = C^~^\x,Cu,fJ')C^'^\y,Cv,l^)'H{x,y,Q,^), 



(27) 



(28) 



at one-loop level. At higher orders, the relation can be a convolution with the functions C's. We 
find then 



T-<{x,y,Cu,Cv,Q,fJ-) 



1 + 



asiK 



I) 



1 

X 



X o_ y O — /^ 1 — ^ 1 —\ 

— m X -\ — In y (In x + In y) 



1 



A* 



A^ 



-I- - In X In -2 + - In y In -2 -I- In -2 -I- In 2 

y 



20Ls , _ , _ 

H m X In y H 

Svr 47r 



„2 g 2 ■ 

/3oln^ - -(3 + lnx + lni/) In ^ 



(29) 



where we only give these terms explicitly: the double log terms and the scale-dependent terms. For 
X — > or y — > we find that in Ji there are no Type I double log's. Therefore, the Type I double 
log's are factorized. 
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4.2 Factorization of Type II double log's 

In this section we will introduce soft factors to factorize Type II double log's. We first discuss 
how to deal the double log from Fig. 2b. As discussed before, in Eq. (13, 14) the term — 1/2 In^ x which 
is divergent with x — > is generated from the region in which the gluon 1 is soft. The contribution 
from the soft gluon can be subtracted by using non-light-like gauge link[19]. The subtraction can 
be done by replacing the quark propagator attached by the gluon 1 with the eikonal propagator: 

1 11 

(30) 



{P2 — — m? + ie —2p2k~ + ie —v-k + ie'' 

with v^^ = (t;"*", w", 0, 0) and 3> f^. Combining factors from numerators one can identify the 
soft contribution as the contribution which is obtained from Fig. 2b by replacing the quark line 
attached with the soft gluon with a gauge link along the direction v. Then one realizes that the 
soft contribution can be subtracted or factorized with the soft factor like 

(0|y(x, -oo, ?;)G"''^(x)G''''^(0)|0), (31) 

V is the gauge link along the direction v from — co to x. At tree- level the soft factor is just the free 
propagator of gluon. This type of subtraction has been discussed in [20J. 

The above soft factor as it stands is not gauge invariant. To factorize the Type II double log in 
Fig. 2b with a gauge invariant soft factor one can modify or extend the above soft factor. For this 
we consider the following soft factor defined with field strength tensor and gauge links: 

S{k) = J d^xe''^ ''Tt{0\T[V{x,oo,w)G-^'{x)V{x,-oo,v) 

V\0, -oo, w)G-' iO)V^{0, oo, w)] |0) {g^, - n^l, - nj^) , (32) 

with k'^ = {xP^ , —K~ ,0,0). The gauge link along the direction v pointing to future or starting 
from the past is defined as: 

V{x,oo,v) = Pexp^—igsJ dXu ■ G{Xv + x)J , 

V{x,—oo,v) = Pexp^—igsJ dXu ■ G{Xv + x)^ . (33) 

The three vectors for directions of gauge links are defined as the following: We take the vector v 
as = {v~^ , v~ ,0, 0). The other two vectors w and w are transverse. The three vectors are chosen 
as: 

V ■ w = 0, V ■ w = 0, w ■ w = 0. (34) 

In the soft factor the limit ^ v~ is taken. The reason for taking the field strength tensor 

is because with other components of G'^'^ the soft factor under the limit 3> v~ can have the 

divergence as 

At the tree-level we have: 

5W(fc) = -z(iV2_i)([^. (35) 

At one-loop level, many diagrams exist. In Fig. 3 and Fig. 4 diagrams for one-loop corrections are 
given. Some diagrams are automatically zero because of the directions of gauge links in Eq.(34). 
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V 



w 




Figure 3: Part of one-loop corrections. 



These diagrams are not drawn in Fig. 3 and Fig. 4. We first study the contribution from Fig. 3b 
which corresponds to the soft contribution in Fig. 2b. The correspondence is reahzed by replacing 
the quark propagator attached by the gluon 1 in Fig. 2b with the gauge link along the direction v. 
It is straightforward to obtain the contribution: 



2 a 2 



(36) 



with 



2v+{K- 



(37) 



From the above result one can see that the contribution from Fig. 3b to the soft factor contains the 
same Type double log as that in the contribution from Fig. 2b to the form factor, as expected. 
The contribution from Fig.3e is exactly zero. The contributions from other diagrams in Fig. 3 

are: 



sm^a 

S{k)\ze 



S{k)\ 



3c 



o + o(C 

0, 



2) 

V I 1 



S{k)Uf = S{k)\sg = sCk)\3h = S^'Hkf'^'^ 



167r 



2 In 



sCk)\ 



3d 



-5W(fc) 



3asNc 



In 



-k^ 



42 



+ 1 



+ 3 



(38) 



There are also one- loop corrections from the self-energy of gauge links, exchanges of gluon 
between gauge links and the self-energy of gluon. The corresponding diagrams are given in Fig. 4. In 
contrast to the contributions from Fig. 3., which do not contain any soft divergence, the corrections 
from the self-energy of gauge links, exchanges of gluon between gauge links in Fig. 4. contain I.R. 
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{a) 



(6) 



Figure 4: Self-energy-Corrections of gauge links. 



singularities. These contributions can be subtracted by modifying the definition of S{k) in Eq.(32) 
For this purpose we introduce 

1 



Sq{x) = ^Tr(0|T \V{x, oo, w)V{x, -oo, v)V''{0, -oo, w)V\0, oo, w)\ |0), 



and modify S{k) to S{k): 



Sik) 



Sq{x)' 



(39) 



(40) 



With the modification, The contributions from Fig. 4a to Fig.4e are subtracted, i.e., there is no 
contribution from Fig. 4a to Fig.4e to S{k). The remaining contribution from Fig. 4 to S{q) is only 
from Fig.4g. It is 



~s{k)Ug = s^'\k)^ 



iV, - In 



31 \ 2^^ / ^' 



^ + 3 



(41) 



The final result for S{q) at one-loop level reads: 



S{k) = S'^^^k) 



^ ln2^ + 21n^ -i^(31iV. + 8iV,)ln^ 



47r V V3 



vr2 + 



9 



(42) 



With the soft factor S we can factorize the Type II double log term in Fig. 2b. 

Now we turn to the Type II double log In^ x in Fig. 2d. It is rather difficult to construct a soft 
factor to factorize the double log. However, at one-loop level we can build the following soft factor 
which can be used for the factorization at least at one-loop level: 



Soik) 



2il ■ k 



d^xe^'^-'^TiiOlT \v{x, oo, n)G-^(x)y(x, -oo, v)V\0, -oo, 10)0-" (0) 
(l)+W(o, 00, w))] |0) {g^, - n^l, - nj^) . 



(43) 
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id) 



Figure 5: The leading order contributions of the correlation function Sd- 



This correlation function consists of a covariant derivative of one gauge link. It is gauge invariant. 
The leading order contribution comes from the gluon change where the gluon is emitted by the term 
of the covariant derivative of the gauge link. In the limit oi v~ ^ the leading order contributions 
to Sd come from the diagrams given in Fig. 5. The contribution from Fig. 5b will contain the double 
log corresponding to that from Fig. 2d. The contribution from Fig.Sd is zero in the dimensional 
regularization. The contributions from each diagrams in Fig. 5. are: 



SDik)\5a 

SD{k)\5b 

SD{k)\bc 
SD{k)\bd 



( 

0. 





In 




Stt 




as 


In^ 






Q 


asN^ 


In 


( 


Stt 





+ 



13 



+ 2 In 



+ 1 



In 



i + o(C') 



(44) 



From the above result one can see that the double log In^ x in the contribution from Fig. 2d. to the 
form factor is correctly produced with the contribution from Fig. 5b. to the soft factor. To factorize 
all Type II double log term of x, we can now introduce the following soft factor combining above 
results: 



\^iS{k) + SDik 



1 



1) 



In^ — ^ + 2 In 



xQ2\ a, / 6 

-{ 19iVc + 8iVf + — 

C2 i 48^ V ^ iVc 



In 



xQ2 
(45) 



where • • • denote constant terms. This soft factor will be used to re-factorize the double log related 
to g"*", i.e., In^ X with q given in Eq.(lO). 

The Type double log of y in the contributions from Fig.2e and Fig.2f can be handle in a similar 
way by introducing another soft factor. The soft factor can be obtained from through the time- 
reversal transformation, where we replace the transformed vector v with u and the transformed 
momentum k with (P"*", — 0, 0). The vector u is given by = (ti^,n^,0,0) with u~ ^ n"*". 
We denote the obtained soft factor as Su{yQ'^ , Cu, fJ-)- Its one-loop result can be read from as: 



Su{yQ'^,Cu,fJ-) 



1) 



C2 



C2 



as 
487r 



19Nc + 8iV/ + ^ ) In 



NJ yQ' 
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2u-{P+f 



(46) 



With the introduced two soft factors we can now factorize the Type II double log's by writing 
the factorized form of the form factor as: 



Svras dx dy 



9Q^ Jo X y 

■7i{x, y, C,u, (^{i, (^uj Cvi Qi A*)) 



(47) 



with 

T~(-{x, y , (^u, (^y, (^u, (^y, Q, ji) 



1 + 



2a. 



1 



In X -\ — In y + In X In -pr- H — In -pr + In 



y 

1 . 2 



yQ"' 



Sm 



1 



+ In 7/ In ^ + ^ In^ ^2" + 1^ ~^|~ + - In x In ^ + - In y In „ 

Cu Cm Cm ^ y Cm 



+ ln-2 +ln-2 

^M ^M 



83a. _ 
+ ^T-; — In — -r + — — In X In y 



247r Q2 



+- 



47r 



8 



^0^"^ ~ ^(3 + lnx + ln?/) In 



(48) 



where we have taken = 2) and Nf = 3. It is clear that the above expression does not contain 
the divergent double log In^ x or In^ y for x — > or y — > 0, respectively. All these divergent double 
log's are contained in the soft factors and NLCWF's. 

5. Resummation and Numerical Result 

To re-sum those double log's in the original H we express the form factor in LCWF's through 
functions C"s and take all scales at Q^: 

FAQ) « ^ /' ^^H^, Q)Hy, Q)C{x, Q, Q)C{y, Q, Q)SAxQ\ Q, Q)Su{yQ\ Q, Q)'Hr{x, y), 
DQ^ Jo X y 

(49) 

with the simple hard part: 



Hr{x,y) = n{x,y,Q,Q,Q,Q,Q,Q) 
2a 



1 + 



Sir 



In X H — In y + In X + In y + In X In y 
X y 



H 



(50) 



so that 7ir{x,y) does not contain those double log's. 

The resummation can be done as in the following. For C-function we can follow |12j to express 
C{x,Q,Q) in term of C{x, fio/^y^, fio) with /xq determined by: 



(51) 
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It should be noted that for a given large scale Q hence a small as{Q) as{fJ-o) is always smaller than 
as{Q) because the asymptotic freedom of QCD. Hence an pertrubative expansion of the C function 
in asipo) still make sense. We have: 



C{x,Q,Q) = exp<^- — 

6po 



Inx — —xasiQ) Inx — X + 1 
4tt 



Inx 



+ 1 



•C(a:,/io/vx,/io) 



C{x,iJ,o/Vx, Ho) 



1 



Qs(mo) 
Svr 



TT 



- 2Li2(x) - — --lnx + — + 2 



+ 0{al). (52) 



All double log's are now re-summed in the exponential factor. 

For the resummation of the soft factor 5^ we have the following ("-evolution equation: 



dSy{xQ'^,Cv,fJ') 2as 



In 



xQ^ 



1 Sy{xQ^Xv,^j)- 



d In Q Svr 

Using this equation we can express Sy{xQ'^, Q, Q) in term of S^{xQ'^\/^Q, Q): 



(53) 



S,{xQ\Q,Q) 



exp 



Svr 



{\n^x + 2\nx^ Sy{xQ'^\^Q,Q), 



S,{xQ^,V^Q,Q) = l-as{Q){---) + 0{al), 



(54) 



where (• • •) contains only constant terms, i.e., no log's. Doing the same for Su we finally have the 
resummed form factor: 



^^<^(x,g)<^(,,Q)e--(^'«)--(-«), 
9 Jo X y 



S{x,Q) = — 
6po 



Inx — ^xUsiQ) Inx — X -|- 1 
47r 



Inx 



+ 1 + 



Sir 



(ln2x-F21nx) , (55) 



since we work at leading order re-summation, we have to take TCr at tree-level, i.e., Tir = 1- 
The above formula is applicable when one knows LCWF's at the large scale Q. If LCWF's are 
determined in a lower scale /i, one can use evolution equation to express (/)(x, Q) in 0(x, fi). In this 
case we have: 



/C(x,y,Q,/x) 



87ras{Q) dx dy 



9 



X y 



-</>(x,;u)</)(y,^)e-^(^'Q)-^(^'«)+^(-*«''^), 



--^ln^(3 + ln(xy-)) 



(56) 



We will use our resummation formula in Eq.(55) and Eq.(56) to give our numerical results. In 
our formulas the nonperturbative input is the LCWF. The LCWF has the asymptotic form if 
goes to oo: 

4>{x, h) = 6x(l - x)/^ H , (57) 

where • • • stand for terms which are zero in the limit fi ^ oo. The LCWF can be expanded with 
Gegenbauer polynomials [T]. A model for (p has been proposed by truncating the expansion |21j: 



(x,^) = 6Ux{l - x) (l + 02(M)cf '(2x - 1) 



(58) 
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Figure 6: Numerical results for the form factor. Curve-^ and Cure-As are drawn by using LCWF 
in Eq.(57) without and with the resummation. Curve- i? and Cme-Bs are drawn by using LCWF 
in Eq.(58) without and with the resummation. 

where (p2i^J'o) is determined by QCD sum-rule method at hq = lGeV[2T]: 



We will use these two types of LCWF to give our numerical results. We will use Eq.(55) with 
the asymptotic form of (j) to make our numerical predictions. For LCWF given in Eq.(58) we use 
Eq.(56). We take the A-parameter as A = 237MeV. The numerical results are represented in Fig. 6.. 
Our numerical results do not strongly depend on the value of A. There is only a little change if we 
change A from lOOMeV to SOOMeV. 

Our numerical results show that the resummation has significant effects. For the asymptotic 
LCWF the difference between the resummed- and unresummed form factor can be from about 30%. 
The difference with LCWF in Eq.(58) is about 50%. We notice that the resummed form factor 
is in general smaller than the unresummed. In [7j a detailed study for numerical predictions at 
one-loop level has been done. With various models of LCWF's it has been shown that the one-loop 
correction is in general positive. The correction is large and can be at order from 40% to 100%. 
This brings up the question if the perturbaive expansion is reliable. Because our resummed form 
factor becomes smaller, we expect that the situation can be improved, or at least partly improved. 
But this needs a detailed study by including higher order corrections. 

In Fig. 6 we also give the experimental result from [22j . In the case of 7*7r — > 7 we have found 
that the experimental data in the region SGeV^ < Q2 < lOGeV^ can be described fairly well with 
the resummed form factor jT2]. But in the case here, we fail to re-produce the experimental results. 
The similar situation also appears in the study at one-loop level in [7], where the experimental 
results can not be reproduced with various models of LCWF. This problem deserves a further 
study. 



M^^o = IGeV) = 0.44. 



(59) 
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6. Conclusion 

The hard part at one- loop level in the collinear factorization for the form factor in 7*7r ^ tt 
contains double log terms of x and y. At n-loop level it is expected that the large log terms like 
In^"" X or In^" y appear. These double log's are dangerous and can result in that the perturbative 
expansion of the hard part becomes a divergent expansion. A resummation of these terms with a 
simple exponentiation can not be done because it results in divergent results. In this work we have 
studied the resummation of these double log terms. For this purpose we have identified the origin 
of the double log's in detail in the first step. In the second step we have employed the concept of 
factorization. We have re-factorized the form factor so that the hard part does not contain double 
log's. In the re-factorization NLCWF instead of LCWF and soft factors are introduced to capture 
these double log terms. 

Non light-like gauge links are used in NLCWF and soft factors. These links play important role 
in the re-factorization and resummation. Because the gauge links are not light-like, the introduced 
NLCWF and soft factors contain extra scales beside the renormalization scale /i. Using evolutions 
of these scales we are able to resum the double log terms into two exponential factors, one is for 
In^ X, while another is for In^ y. An interesting aspect of our approach for the resummation is that 
the resummed form factor only contains LCWF's as nonperturbative quantities. Every ingredient 
in the resummation except LCWF's can be calculated perturbatively. 

With the resummation we have given numerical results for the form factor with two choices 
of LCWF. We have found that the resummation has significant effects. Between the resummed- 
the unresummed form factor the numerical difference is at level about 30% — 50%. With various 
models of LCWF's it has been shown in |7| that the one-loop correction is in general positive. 
The correction is large and can be at order from 40% to 100%. Because of the large correction, 
the perturbative prediction can be unreliable. With the structure of the exponentials for the 
resummation one can find that the resummed form factor becomes smaller than the unresummed. 
With this fact we expect that the situation with the perturbative expansion can be improved, or 
at least partly improved. But this needs a detailed study by including higher order corrections. 

In the case of 7*7r ^ 7 we have found that the experimental data in the region 3GeV^ < < 
lOGeV^ can be described fairly well with the resummed form factor [12]. It is frustrated in the 
case of 7*7r — > vr if one compares theoretical results of collinear factorization with experimental 
data. With or without resummation, the experimental data at higher is not in agreement with 
perturbative results. This deserves a further study in experiment and in theory. In this work 
we have performed the resummation at one-loop level. It is possible to extend our work to the 
resummation of the remaining single log terms and beyond one-loop level. 

Appendix: Gauge Variance in the /c^^-Factorization 

In the collinear factorization the transverse momenta of partons entering hard scattering are 
neglected at leading twist. As we have seen that the perturbative part in Eq.(9) contains double 
log's at higher order and is divergent when the momentum fraction x or y approaches zero. It has 
been suggested that for small longitudinal momenta one may need to take the transverse momenta 
into account. This leads to the so-called fcT-factorization. The /c^-factorization has been widely 
used in exclusive i?-decays. Although it has been widely used, the factorization has not been studied 
beyond leading order except the case with tt'j* — > ^|23j in which the study is perform with Feynman 



18 



gauge. Because the perturbative parts in the factorization are extracted from scattering amphtudes 
of off-sheh partons and scattering amphtudes of off-sheh partons are not gauge-invariant, it is not 
expected that the perturbative parts are gauge- invariant. In [15] it has been pointed out that such 
a factorization is gauge-dependent because the perturbative parts contain soft divergences at loop- 
level which depend on gauges. It should be noted that in the case with tt^* — > 7 the problem of 
gauge invariance appears beyond tree-level, because gluons are exchanged at one- or higher loops. 
Taking the case with tt^* — > vr, it is easy to show the problem at tree-level, because at that level 
gluon-exchange happens. 

In the A;T-factorization, the form factor in Eq.(2) can be factorized in a similar way as in Eq.(9): 

FiriQ) ^ dx(fki±dy(fk2±(l){x,ki_L)(j){y,k2_L)H{x,y,ki_L,k2±). (60) 



The definition of the /cT-dependent wave functions can be found in [151 HTJ [23] . The perturbative 
part H is extracted from scattering of off-shell partons. Instead of parton momenta in Eq.(4) one 
has the momenta for these off-shell partons: 

p1 = {xoP^,0,P±), = (xoP+,0,-px), ki^ = (o,yoK-,kj_) , kl^ = (o,yoK- , -k^) . (61) 

At leading order the same diagrams in Fig.l. give contributions to the form factor and hence to H. 
In calculating these one makes the projection by replacing the spinor products like those in Eq.(5) 
with: 

u{pi)v{p2) 7^75, v{ki)u{k2) 1^1^ ■ (62) 

We take /i = — as before. Then the contribution to the form factor from the n-quark is only from 
Fig. la.. If the /cr-factorization respects the gauge invariance, this contribution must be gauge- 
invariant. In Feynman gauge one finds the hard part from the n-quark: 

H{x,y,p^,ki_)\u = —- ^ (63) 

xyQ^ + {p±- k±y 

In deriving this result one has used the power counting: xQ ~ yQ ^ p± ^ k± ^ 5 and only the 
leading term in 6 has been taken into account. To extract H one has also used the tree- level result 
of fcr-dependent wave functions with the off-shell partons: 

4'{x,ki±) 6{x - xo)6^{ki± -pj_), 4'{y,k2±) 6{y -yo)6'^{k2±-kj_). (64) 

It should be noted that there is no problem of gauge-invariance for the tree-level result of wave 
functions, because there is no gluon exchange. 

Now we calculate the hard part in an axial gauge fixed with an arbitrary vector w, i.e., w-G = 0. 
The gluon propagator in this gauge reads: 



7^ -|- ie 



g'^ h w 



w ■ q {w ■ q) 



(65) 



In the above the first term is just the propagator in Feynman gauge. We obtain the contribution 
from the n-quark to the hard part in the axial gauge as: 



TTf h \\ - I Ji 'w±^{p-k)± 2 P±-iP- k)± 

H[x,y,p^,k^)\u - --Q2^^j^^_j:^yV w.{p,-k,)^'' 2{w.{p,-k,)f 

' ^^-^ (66) 



2n; • (pi - ki) 



—r^k±_ ■ [k - p±) + w ■ k_L 

yK 
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In this result the momentum fraction xq in pi and yo in in Eq. (61) should be replaced with x 
and y, respectively. In the above {• • •} every term is at order of 0{d^). Therefore no term can be 
neglected with the power counting. This result clearly indicates that H is gauge-dependent. It is 
interesting to note that the first two terms are canceled and all terms in H depend on w, if we take 
the vector w as = (0,0, u;^). For showing this clearly we give the explicit result for this case: 

H{x,y,p±,k±)\u = —— — — ^^{tTi 7 Tr\y2+^ 1 TT f ' 



One can also show that the hard part is gauge-dependent in the general covariant gauge [24]. In 
the gauge the gluon propagator reads: 



4- ie 



- 



In this gauge the contribution from the u-quark reads: 

H{x,y,p±,k±)\u = —- — 1 + — ^ • (69) 

xyQ^ + (p± - k^y V 2 xyQ^ + {p± - ki_Y J 

Here the gauge-dependent term in (• • •) is proportional to a. This term is at order of 0{6^) 
in comparison with the first term. Therefore, the gauge-dependent term can not be neglected. 
One notices here that the transverse momenta appear in the numerator. One may argue that these 
terms with the transverse momenta in the numerator may be factorized with higher-twist operators 
other than the leading- twist operator used to defined 4'ixyPi.)- Even if one can do so, these terms 
contributing to the factorization hence to the form factor are still gauge-dependent and can not 
be neglected with the power counting in comparison with the term factorized with (p(x,p±). The 
conclusion here is that the /cr-factorization here at tree- level is gauge-dependent. Beyond tree-level, 
the hard part will receive gauge-dependent contributions which are divergent |15i I16j . 
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